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TECHNICAL NOTE 3105 

AERODYNAMICS OF SLENDER WINGS AND WING-BODY COMBINATIONS 
HAVING SWEPT TRAILING EDGES 
By Harold Mirels 


SUMMARY 

A general method, based on two-dimensional crossflow concepts, is 
presented for obtaining the lift and moments on highly swept wings. 
Emphasis is placed on obta ining solutions for wings having swept 
trailing edges. The method is applicable for all problems where the 
velocity boundary conditions can be made homogeneous by differentia- 
tion in the streamwise or spanwise directions. 

Lift, roll, and pitch solutions, for highly swept wings, are pre- 
sented. Both direct problems (where the plan form is given) and in- 
verse problems (where the shed vortex sheet is given) are considered. 
The solutions of the direct problems are expressed in terms of func- 
tions which are evaluated from integral equations. Some limiting 
solutions of the integral equations are indicated. Numerical results 
are given for wings having parallel leading and trailing edges. 

The transformation of a wing-body problem to an equivalent iso- 
lated wing problem is discussed and the application for finding the 
lift of a wing-body combination is indicated. 

Application of the method for • solving unsteady two-dimensional 
inconpressible flow problems is also indicated. In particular, the 
Wagner problem is formulated in terms of the techniques developed 
herein. 


INTRODUCTION 

In 1924, Munk (ref. l) published a remarkable paper concerning 
the calculation of the aerodynamic forces on airships. His theory 
was based on the idea that the velocity field induced by a slender 
body is essentially two-d im ensional in planes transverse to the 
body axis . This reduced the complicated three-dimensional problem 
to an equivalent two-dimensional unsteady-flow problem and permitted 
the use of very elegant methods of solution. 
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With the advent of transonic and supersonic flight, interest in 
slender bodies was renewed. In 1946, Jones (ref. 2) revived Munk's 
ideas and used them to compute the forces on low-aspect-ratio pointed 
wings. He indicated that for such wings, compressibility has no ef- 
fect and his results applied equally at both subsonic and supersonic 
speeds. Following Jones' ■ example, many papers were written on the 
aerodynamics of slender wings, bodies, and wing-body combinations. 

See, for example, references 3 to 11. One of the most notable of 
these is Ward's paper (ref. 5) which provided a rigorous justifica- 
tion for the Murik- Jones approach. Ward showed that their solution 
may be considered as the first term of an expansion in terms of a 
"slenderness" parameter. Ward considered supersonic flight speeds. 

More recently, Adams and Sears (ref. 10) gave a similar result for 
the subsonic case. 

Relatively few investigations have been made for cases where the 
two-dimensional crossflow contains a shed vortex sheet. This occurs, 
for example, when a slender wing has a swept trailing edge. In this 
case the crossflow generally contains a shed vortex sheet of unknown 
strength. The crossflow is not independent of upstream conditions 
and the problem is considerably more complicated than those considered 
by the early investigators. References 6 to 8 considered wings having 
swept trailing edges. In these references the distribution of vorticity 
in the shed vortex sheet is assumed and the wing plan form which would 
give rise to such a distribution is then found. This is the so-called 
inverse problem of aerodynamics and is generally simpler to solve than 
the direct problem. The direct problem is one in which the wing is 
completely specified and the flow field is to be determined. Robinson 
(ref. ll) appears to be the only one to have presented a solution of 
the direct problem for wings with swept trailing edges- 1 -. He has 
treated the lift problem. In reference 12 it is indicated that 
Robinson's solution is applicable only when the trailing edge is 
slightly swept. 

In the present paper, a general method is developed for solving 
low-aspect-ratio problems involving shed vortex sheets. Both direct 
and inverse problems are considered. The method is applicable for 
all planar problems where the velocity boundary conditions can be 
made homogeneous by differentiation in either the streamwlse or the 


-*-A private communication frcm K. W. Mangier, in connection with 
reference 12, served to call the author's attention to reference 13. 
In reference 13, Mangier has independently obtained the lift, roll, 
and pitch solutions for a highly swept wing (direct problem) by a 
method similar to that used herein. A comparison of reference 13 
with the present report is given in appendix E. 
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spanwise directions. The lift, roll, and pitch solutions, for a highly 
swept wing, are presented and some numerical examples are worked out. 
Applications to wing-body interference problems and to unsteady two- 
dimensional incompressible flows are also indicated. 

The research reported herein was conducted at the Graduate School 
of Aeronautical Engineering, Cornell University. Thi s paper is based 
on material which was originally presented to the faculty of the 
Graduate School of Cornell University, in June 1953, as a thesis for 
the degree of Doctor of Philosophy. The author wishes to express his 
sincere gratitude to Professor N. Rott for his advice and criticism 
during the course of the study. The author also wishes to thank the 
other members of the staff of the Graduate School of Aeronautical 
Engineering and Drs. H. K. Cheng and M. C. Adams (both formerly at 
Cornell University) for stimulating discussions. 


SECTION I - BASIC CONCEPTS 


In the following section, the basic results of previous investi- 
gators, primarily reference 5, are summarized. 


The assumption of slenderness is introduced into the equations of 
motion. The general features of the crossflow and formulas for lift, 
drag, and moments are discussed. Finally, symmetry of the velocity 
components, in planar problems, is mentioned. 


1.1 Equation of motion . - Consider a slender body in a free stream 
of velocity Uq, Mach number Mq, pressure Pq, and density p Q (fig. l) . 

The coordinate system is stationary with respect to the body and is de- 
fined such that the x-axis is parallel to the free stream (i.e., wind 
axes) . 


If the body is assumed to perturb the main stream only slightly, 
the equation of motion can be linearized and reduced to the Prandtl- 
Glauert equation: 



d 2 (p 3jP = 0 

cbc^ dy 2 Bz 2 


( 1 . 1 . 1 ) 


where (p is the perturbation velocity potential. Thus, if u, v, and 
w represent the perturbation velocities in the x, y, and z direc- 
tions, respectively, then 
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Equation (l.l.l) applies for both subsonic and supersonic flight. To 
as good an approximation, the pressure at any point in the flow field 
is given by 



r Sep 1 / 

a<A 2 (<■ 

><A 2 ‘ 


P - Pq Pq * 

^0 ee + 2 || 

§ y) + V 

sjj 



( 1 . 1 . 2 ) 


where Pq and p 0 are the pressure and density, respectively, in the 
undisturbed flow. 


A body is considered slender if 



(1.1.3) 


where c 0 characterizes the length of the body and bg characterizes 
its width. Under this condition, equation (l.l.l) becomes 



Sy2 ^ z 2 


(1.1.4) 


which is the governing equation in slender body theory. Equation 
(1.1.4) indicates that 9 can be found, to within a function of 
x, by considering the flow in each yz-plane to be a two-dimensional 
incompressible flow. 

The boundary conditions for equation (1.1.4) are usually expressed 
in terms of the perturbation velocities. Let v and p be orthogonal 
coordinates which are normal and tangential, respectively, to the curve 
defining the body cross-sectional area in a particular yz-plane (fig. 

2) . Let v n = dcp/dv be the perturbation velocity in the v-direction. 
Then, the condition that the resultant flow be tangent to the body sur- 
face requires that, at the body surface (ref. 5), 

v n = U 0 ^ ' (1.1.5) 


1.2 Asymptotic form of crossflow . - Two-dimensional incompressible 
flows are well understood and are most easily handled in terms of the 
complex variable £ = y + iz. Let W = cp + it where W is the com- 
plex potential function, <p is the velocity potential, and t{r is the 
stream function for the crossflow. (The functions W, <P, and con- 
tain x as a parameter.) The Laurent expansion for W, valid every- 
where outside the smallest circle, center at £ = 0 , enclosing all 
the singularities of the flow, is 
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¥ = f + f o In 



m=l 



( 1 . 2 . 1 ) 


where f, ±q, f^, . • • axe functions of x. The function f intro- 
duces a constant pressure at each crossflow. Formulas for its evalu- 
ation are given in reference 5 (for Mq > l) and reference 10 (for 

Mq < l). Since f does not contribute to lift, and since only the 


lift problem is of interest here, those formulas will not be repeated. 
The function fQ is generally complex. The real part is proportional 


to the source strength required to simulate the expansion or eontrac- 

Uq . dA cs 

tion of the body and equals — (ref. 5) where A cs is the cross- 

r 


2 it dx " cs 

sectional area of the body. The imaginary part equals -r/'2it where 
is the net circulation at each section. Equation (1.2.1) can then 
be written 


¥ 




^cs 

dx 



m=l 


( 1 . 2 . 2 ) 


The complex velocity V = v - iw is found by differentiating 
equation (1.2.2) with respect to £ and equals 


•o 

v * s ( u ° 75T - ^ f fin E ' (m+1) (1 ’ 2 ' 5) 


m=l 


Differentiation of equation (1.2.3) with respect to £ gives 


bV 

5? 


h ( u ° 


2f] 


6f 2 

t 4 


(1.2.4) 


Similarly, if equation (1.2.3) is differentiated with respect to x 


bV 1 ( d2A cs \ i df x /dx df 2 /dx 

s ’ 2? V 0 ~^r) X ' ~lz~~ ' “F - ' ' ' ' 


where it has been assumed that dr/dx =0 so as to satisfy the law of 
conservation of circulation. Equations (1.2.2) to (1.2.5) define the 
asymptotic behavior of ¥ and its derivatives. These expressions 
will be useful in later developments. 



6 


NACA TN 3105 


1.3 Lift, drag, and moments . - Expressions for computing the forces 
acting on a body are summarized herein. They are based primarily on 
the results of reference 5. It will be assumed that r = 0, which is 
the case for al 1 problems treated in the present report '(except for 
section 2.3). 

The net force which has acted on a body upstream of any section 
x = constant is obtained by a contour integration about the body 
(fig. 2) and equals 


F + 

y 


IF, 



<P cl£ 


(1.3.1) 


Equation (l.3.l) can be evaluated by Cauchy's theorem for the path 
Cf + Cg indicated in figure 2. The cut is introduced to make W 

single-valued. The path C 2 is sufficiently far from the body so 

that the Laurent expansion for W can be used along this path. 
Since there are no singularities inside the path + C 2 , Cauchy's 
theorem gives 
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where is the point common to paths Cp and C2, as indicated in 

figure 2. Let £g = yg + i z g be the center of area of the cross- 
sectional area A cs . Then, by definition of the center of area 

£ dy dz 


Let v and p be coordinates normal and tangential to the path Cp 
and let v n be the local velocity component normal to the path. Then 


' A = 
>g cs - 



_d_ 

dx 


^g -^cs) 




d y , 

? S ^ 


_1 

U 0 


Thus, equation (1.3.2) can be expressed in 




£ a* 


f l 


U, 


0 d 


F y + - 2 «eo u o |'i*SS < A cs 5gJ 


The force per unit x is 
dF, r dF 


d5T + 1 d5T = 2jtp 0 U 0 


fal 

dx 


u 0 d 2 


+ 2jt { = bc 2 ^ A cs £g) 


(1.3.3) 


(1.3.4) 


The net moment about the z- and y-axes is then given by 


M Z + iMy = 




(1.3.5) 


which can be evaluated in terms of equation (1.3.4). Equations (1.3.3) 
to (1.3.5) were first derived in reference 5. 
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A similar expression can be obtained for the rolling moment about 
the x-axis. Frcm consideration of impulse, the net rolling moment 
which has acted on a body up to a certain section x is 



<p(y <3y + z dz) 


(1.3.6) 


The author was not able to evaluate equation (1.3.6) by a procedure 
similar to that used for equation (1.3.1) except for the case of zero 
thickness. Thus, when dz = 0 in the integrand of equation (1.3.6), 





r \ 

M x = R.P. 

-poUo © 

= R.P. 

-POUo © (wd( 


Jc 1 


Jc 2 


= I.P. (2*p o U0f 2 ) (1.3.7) 

where R.P. and I.P. indicate the real and imaginary parts, respec- 
tively. The rolling moment, per unit x, for this case, is then 

sr - I - p - ( 2 ”«>oo sr) f 1 - 3 - 8 ) 

From equations (1.3.3) to (1.3.8) it is seen that the forces and 
moments acting on the body can be obtained from a knowledge of the 
asymptotic form of W and the cross-section geometry. The results 
are equivalent to integrating the pressure distribution over the 
body surface, providing the quadratic terms in equation {1.1.2) ape 
retained. The use of equations (1.3.3) to (1.3.8) is considerably 
simpler them the surface integration of pressures. 

For lifting surfaces, there is an induced drag associated 

with the lift, which can be determined from a consideration of the 
suction force along the airfoil leading edge. Thus 


Di 



dF„ 


(1.3.9) 


where F x is the suction force acting in the positive x-direction 
and the contour integral is taken along the airfoil leading edge. 
(The symbols L and F z will be used interchangeably to indicate 

net lift.) The magnitude of the suction force is obtained by con- 
sidering the flow at a subsonic leading edge to act locally like 
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that around a two-dimensional airfoil (fig. 3) . In figure 3, the sub- 
scripts n and t represent components normal anri tangent to the 
leading edge. The suction force per unit span on the leading edge 
of a two-dimensional airfoil in a free stream of Mach number Mq n 
is ' 


dx t 


= - P 0 rt 




“ 5,1 


Tim 
x n*0 




Assuming this to apply loca l ly, for a swept wing in a free stream of 
Mach number Mq, 

= " P 0 Wl " MO 2 sin 2 0 lim % 2x n 

x n ->0 

But, from figure 3, 


«g n _ ^x 
dxt dy 2 


v 

Un cos 0 

Xn = &2 ~ y)cos 9 


Then 


££x 

dy 2 


pq* 



- Mq 2 sin 2 0 
cos 9 


lim 

J-*y 2 


T 2 (y 2 - y) 


(1.3.10) 


Up to this point, the slender body assumption has not been made in 
equation (l.3.10). Since v in equation (1.3.10) will be obtained 
from slender body theory, it is consistent to expand 

vr - Mq 2 sin^0/cos 0 in terms of the slenderness parameter. Thus, 




Mq 2 sin 2 0 
cos 0 


= 1 + 0 


Ia/mq 2 - jj 


c 0 


where 0( ) indicates order of magnitude. For a slender configuration, 
equation (1.3.10) becomes 
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^2 


= - p 0 ot lim v 2 (y 2 - y) = - 

y -^2 


pp* 

(dy 2 /dx) 2 


11m 

y^2 


U 2 (y 2 - y) 

(1.3.11) 


The suction force term in equation (1.3.9) is obtained by integrating 
equation (l.3.1l) along the airfoil leading edge. An alternate method 
for computing Dj_ is to find the kinetic energy in the wake behind the 

wing. Both methods give the same result. 


1.4 Symmetry in planar problems . - The solution of aerodynamic 
problems is greatly simplified by imposing, at an early stage, what- 
ever symmetry the velocity field must have. The symmetry which ex- 
ists in problems involving zero thickness wings will now be noted. 


If a wing has moderate camber and twist and is at a small angle 
of attack to the main stream, the solution can be obtained by speci- 
fying the boundary conditions in the z = 0 plane rather than on 
the surface of the wing. The w velocity is symmetric while the 
u and v velocities are antisymmetric about the z = 0 plane. 

That is. 


u(x,y,z) = - u(x,y,-z) 
v(x,y,z) = - v(x,y,-z) } 
v(x,y,z) = w(x,y,-z) 


(1.4.1) 


Thus the w velocity is continuous while the u and v velocities 
are discontinuous (equal and opposite) or zero for corresponding 
points on the upper and lower surface of the z = 0 plane. 


When the w boundary condition on the wing is symmetric with 
respect to the y = 0 plane (as in the case of a lifting or pitching 
wing) and the corresponding plan form edges have the same type singu- 
larities, then the following symmetry also applies: 


u(x,y,z) = u(x, -y, z) 
v(x,y,z) = - v(x, -y, z) > 


w(x,y,z) = w(x, -y, z) 


(1.4.2) 


When the w boundary condition on the wing is antisymmetric with re- 
spect to the y = 0 plane (as in the case of a rolling wing), then 
equations (1.4.2) are replaced by 


3120 



CY-2 tack 


MCA TN 3105 


11 


cn 

H 

DO 

o 


u(x,y,z) = - u(x,-y,zj 

v(x,y,z) = v(x,-y,z) \ (1.4.3) 

w(x,y,z) = - w(x,-y,z) 

again assuming that corresponding edges have the same s ingular ities. 

The symmetry of the derivatives of the perturbation velocities is 
found by differentiating equations (1.4. l) to (l.4.3). 


SECTION 2 - GENERATING FUNCTIONS 

The functions and Sv/Sx axe termed generating functions 

herein since they are not of particular interest in themselves, but 
their integration leads to the solution of flow problems^. Generating 
functions can be used to solve al 1 planar problems for which the bound- 
ary conditions in the z = 0 plane are made homogeneous by differenti- 
ation in the x- or y-direction. Expressions for the generating func- 
tions axe derived in the following section. 

2.1 Evaluation of branch points . - The functions which arise in 
wing problems usually have branch points on the y-axis. It is essen- 
tial to develop a systematic procedure for introducing cuts so as to 
make these functions single-valued. In all cases, the cuts will be 
introduced along the y-axis to the left of the branch points. 

As an example, consider the function (£ - y n )^ + -'-/2 where N Is 
any integer (0, +1, +2, . . .) and y n is an arbitrary point on the 
y-axis. This function has a branch point at y n . To make the func- 
tion single -valued, a cut is introduced along the y-axis from - -» 
to y n (fig. 4) . Define the function to be real and positive for 


^The application of generating functions to aerodynamic problems 
was first brought to the author's attention by Dr. H. K. Cheng who 
used a similar approach in his thesis "Thin Wings in Conical Flow," 
Cornell University, 1952. The use of the generating function dv/d£ 
is a classical approach to solv ing Laplace's equation, for a certain 
class of boundary conditions, by means of the complex variable. To 
the author's knowledge, the application of the generating function 
dv/dx does not appear explicitly in the literature except for 
Mangier' s use (ref. 13) of an equivalent function du/d£, (see 
appendix E) . 
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£ = y > y n . To evaluate it along the cut, let £ = y n + ee*® and let 
6 go from 0 to ±it. The result is 

(£ - y n ) N+l/2 = (y - y n ) H+1 / 2 for £ » y > y n 

= ±i(y - y n ) N Vy n - y for £ = y< y n ( 2 . 1 . 1 ) 

Thus the function goes from a purely real function for £ = y > y n to 
a purely imaginary function (discontinuous across the cut) for 

5 = y < y n . 

2.2 Behavior of flow near boundary edges . - A boundary edge is 
defined herein as a point in the z = 0 plane where the boundary 
condition changes from a specification of v to a specification of 
w. Such a point generally corresponds to the edge of a wing panel. 

Consider a boundary edge y n with the boundary conditions 
w = 0 for y > y n and v = 0 for y < y n (sketch l). 


z 

a 

v = 0 


w = 



0 


♦►y 


Sketch 1 


These boundary conditions are satisfied by V = (£ - y n )^ + -*-/^ where N 

is an integer. (Because of the uniform procedure prescribed for intro- 
ducing cuts, this solution for V should be coiisidered as valid for 
the upper hkli plane z 0. The solution for the lower half plane 
is obtained from symmetry considerations . ) For physical reasons the 
smallest permissible value of N is -1. A smaller value would make 
W infinite at y n and this is unrealistic. If the v and w 

velocities are to be continuous at y n (Kutta condition), the minimu m 
permissible value of If is 0. 
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The boundary conditions can now be generalized so that w = w(y) 
for y > y n and v = v(y) for y < y n (sketch 2). In the immediate 


z 

* 


W = v(y) w = w(y) 


■►y 


yn 


Sketch 2 

vicinity of the edge, V must have the form 

V ~[j(y n ) - iv(y n )J + B(£ - y n ) n+1 / 2 (2.2.1) 

where B is a real constant. The first term of equation (2.2.1) sat- 
isfies the inhomogeneous boundary condition at y n while the second 

term satisfies the homogeneous boundary conditions indicated in sketch 
1. The exponent N has the same limitations as for the flow in 
sketch 1. 

Thus, the homogeneous part of V has at most a half- order singu- 
larity at an edge and behaves like 

V ~ 7 _ 

<£ - y n ) 1/s 


^ ‘ ' 2 « - y n ) 5 / 2 

where B is real or imaginary, depending on whether the w boundary 
condition is specified for y > y n or y < yn- If the Kutta condi- 
tion applies, the homogeneous part of V behaves like 

V ~ B(£ - yn ) 1 / 2 

dV 1 B 

2 (£ - y ^) 1 / 2 


(2.2.3) 
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Differentiation of V with respect to t, increases the order of the 
s ing ularity at y n . If y n and B are considered functions of x, 

differentiation of equations (2.2.2) and (2.2.3) yields, respectively, 


* +i »\&) 
s 2 (E - y n ) 3 / 2 


Sv 

SE 


1 

2 


B 

(T 



(Kutta) 


(2.2.4) 


where only the leading term is retained in each expression. Thus, 
differentiation with respect to x also increases the order of the 
singularity. If y n does not vary with x, the derivatives in 

equation (2.2.4) become, respectively. 


dB 

3V dx 

(£ ~ Yn) 1 / 2 


§E " S « - Jn) l/2 CK^tta) 

so that the order of the singularity is not increased by the 
differentiation . 


2~5 Determination of generating functions . - 

(a) Special class of flows : In order to lead smoothly to the 

discussion of generating functions, it Is convenient to solve the 
two-dimensional incompressible flow associated with one or more 
flat plates in uniform translation. 

Consider a flat plate to be, at a given instant, on the y-axis 
between y^_ and yg and to be moving vertically downward with the 

velocity w = - oUq (fig. 5) . The boundary con d itions in the z = 0 
plane are indicated in the figure. The complex velocity must equal 


1 ±A 


ys - Vi vr 


+ iaUo 


*2, 


V = - ixxUo 


(2.3.1) 
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where A is an arbitrary real constant which defines the net circula- 
tion. (From the asymptotic form of equation (2.3.1) it can be shown 
that T = 2mU 0 [A + (y^ + y 2 )/2].) 

Equation (2.3.1) is a well-known result. It can be constructed, 
by inspection, in the following way. The second term on the right- 
hand side satisfies the nonhamogeneous boundary condition w = - aUQ 

for yq < y < y 2 • The bracketed term on the right-hand side must 

therefore satisfy the homogeneous boundary conditions w = 0 for 
y^ < y < y 2 and v = 0 for y < y^_ and y > y 2 . From equation 

(2.2.2) it is known that V can have half-order singularities at 
y^ and y 2 . These are introduced as the product i/<>/£ - y^ - y 2 
so that, in the z = 0 plane, the term is purely real for yj_ < y < y 2 
and purely imaginary for y < yq and y > y 2 , thereby satisfying the 
required homogeneous boundary conditions. The factor <xUq(£ + A) is 

then introduced into the numerator to satisfy the condition that V 
behave like l/£ for £ -*■ -. 

The permissibility of introducing A into equation (2.3.1) cor- 
responds to the fact that, mathematically, the circulation about a 
given airfoil is arbitrary. Taking A equal to -y^ or -y 2 is 

equivalent to applying the Kutta con di tion at or y 2 , respec- 

tively. Either choice gives the same value of r, but with differ- 
ent signs. The net circulation is zero when A - - (y^ + y 2 )/2. 


Equation (2.3.1) can. be generalized. This is the problem of m 
wing panels, each moving downward with velocity w = - aU Q (fig. 6(a)). 

There are 2m boundary edges and the solution is 



where the A n are real. The construction of equation (2.3.2) is simi- 
lar to that of equation (2.3.1). The second term on the right-hand 
side satisfies the nonhamogeneous boundary condition w = - cxUq for 

points on the wing panels. The bracketed term satisfies homogeneous 
boundary conditions in the z = 0 plane, (w = 0 on the wing panels 
and v = 0 off the wing panels) since it is purely real for points 
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on the -wing panels and purely Im aginary for points off the wing panels. 
The numerator is a polynomial of order m since this is the highest 
order polynomial which will still satisfy the boundary condition 
V ~ l/i for i -+ -. 


Equation (2.3.2) has m arbitrary constants, which corresponds 
to the fact that the circulation about each of the m wing panels is 
arbitrary. The constants are determined by (l) specifying the circu- 
lation about each wing panel, (2) specifying the Kutta condition at 
m edges, or (3) combinations of (l) and (2). If the Kutta condition 
is specified at m' edges, designated by y^, and the remaining m" 

edges are designated by y^, then equation (2.3.2) can be written 


V = 


iaU 0 i 


m' 


(m-m' -l) 

TT -Jt-K 


^ (m-m ’ ) + A ^n 

Ln=l J 


L n =0 J 


m" 

TT VT^ 


n=l 


> + iaU 0 

(2.3.3) 


Note that m 1 < m and m 1 + m" = 2m. There are now (m - m ! ) constants 
to be determined by satisfying the circulation boundary conditions. 

(if the problem has symmetry in respect to y = 0, the number of un- 
known constants can be reduced by inspection. ) 


(b) The generating function 5v/5£ : The boundary conditions in 

the z = 0 plane can be further generalized to the case where the 
v and w velocities are constants on segments of the y-axis (fig. 
6(b)). The problem is one with homogeneous boundary conditions, both 
in the z = 0 plane and for £ ~y •*, if rather than V is 

considered. The solution for can be found by inspection 

and equals 


dV 



(2.3.4) 


where the A n are real. The denominator of equation (2.3.4) has 
3/2-order singularities since differentiation of V with respect 
to £ increases the order of the singularity at an edge (eq. (2.2 2)). 
The boundary conditions <dv/c) y = 0 and 3w/c3y = 0, alternately, in 
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i 


to 

i 

fcH 

o 


the z = 0 plane, are satisfied, since equation (2.3.4) becomes imagi- 
nary and real, alternately. The boundary condition at infinity is sat- 
isfied since the leading term of the asymptotic expansion of equation 
(2.3.4) is of the form l/£^, which is appropriate for problems having 
a net circulation. There are 3m - 1 constants in equation (2.3.4) of 
which only m are arbitrary. The 3m - 1 constants correspond to the 
m specified w boundary conditions on the wing panels, the m - 1 
specified v boundary conditions in the segments between the wing 
panels, and the m arbitrary circulations which can be imposed. The 
circulation boundary conditions can be replaced by the Kutta condition 
(now a half-order singularity) at m, or less, edges. If the Kutta 
condition is imposed at m* edges (designated by y^) and the remaining 

m" edges are designated by y^, equation (2.3.4) becomes 


3v 

^ = 


(3m-m' -2) 

Z 


n=0 


Tt « - yi) l/2 l Ft cc - yy 3/2 

n=l - ’ 


_1 L- 


m" 

Z 


(2.3.5) 


As a further extension, the wing panels may have constant downwash 
specified but with a finite number of jumps in w occurring across 
each wing panel (see fig. 6(c)). The solution can be found, from the 
previous equations. If each jump is artificially separated by a vortex 
sheet (with the Kutta condition applied at each edge of the sheet) and 
the limit Is then taken as the intervening vortex sheet width goes to 
zero. Similarly, if there is a discontinuity In v at a point in a 
vortex sheet, an intervening wing panel is introduced and then made 
to go to zero. Thus, if such discontinuities occur at m m points 
(designated by y ”* ) , equation (2.3.5) becomes 


l/2 (m J +3m"+2m"' -4) 


Z 


_ J 

~ 1 

n=0 


ft (s - yA ) l/2 

TT (S - y'n ^ 2 

1 

it= 

i 


n=l 

n=l 

L n=l J 


> 


(2.3.6) 


It is noted that the discontinuities introduce singularities of order 
one in the expression for <dV"/5£- This is due to the fact that they 
correspond to log(£ - y” 1 ) type terms in the expression for V. 
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A final generalization will be introduced. Equation (2.3.6) sat- 
isfies homogeneous boundary conditions, for c*V/d£, both in the z = 0 
plane and at «•. A nohhamogeneous constant boundary condition in the 
z = 0 plane can be readily incorporated. Thus, if the w boundary 

'bw\ 


condition on the n^ wing panel is w = C n + 


y, where 


0 


(Do 


is 


a constant and the same for al 1 wing panels, equation (2.3.6) becomes 


SV 


-i< 


l/2(m , +3m"+2m ,n -2) 

^l/2(m , +3m"+2m ,, ‘) + A ^n 


n=0 


r TT (t - yi) 1/2 lf ft « - yS) 3/2 Tft tt - y ff) 

n=l 


n=l 


r“™n« 


n=l 



>- i 

(2.3.7) 


where the additional constant introduced into the numerator of the 
bracketed term is e limina ted by setting equal to zero the coefficient 
of the l/£ term in the asymptotic expansion. 


The inverse problem of slender wing theory is to find the trailing 
edge corresponding to a given shed vortex distribution. The generating 
function dv/d£ will be used, in later sections, to solve such 
problems . 

(c) The generating function dv/chc ; .In the previous paragraphs 
the v or w velocity was specified for all points in the z = 0 
plane. However, in the direct problem of slender wing theory, the 
w is specified on the wing but the v distribution off the wing 
is generally unknown. Consider the case of a swept wing (fig. 7) 
at angle of attack a in a free stream of velocity Uq. Behind 

the trailing edge there is a shed vortex sheet whose strength varies 
only with y. The crossflow then contains two wing panels having the 
boundary condition w = - oUo with an intermediate vortex sheet of 
unknown strength. Previously V was differentiated with respect 
to £ in order to obtain a homogeneous boundary- value problem. It 
is apparent that, for the direct problem, differentiation with re- 
spect to x will yield the same result. Thus dv/dx has homogene- 
ous boundary conditions and the solution for Sv/5x can be con- 
structed in terms of the singularities at the boundary edges. 


Differentiation of V with respect to x increases the order 
of the singularity at an edge, provided the edge varies with x 
(section 2.2). Thus, for a multiwing panel problem having the 
boundary condition w = f(y) + g(x) specified for the wing panels, 
and v = v(y) off the wing panels, the generating function is 
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dV 

EE = 


i< 


dg(x) fl/2(m'+3m"+2m'») 
dx b 


l/2(m'+3m"+2rn m -2) 

+ Z 


A n C n 


n=0 


TT C£ - 

IT (5 - y£) 3y/2 

TT (E - y£’) 

n=l 

n=l 

n=l 


>- 


4 ^g(^) 

dx 

(2.3.8) 


where the notation is the same as that used in equation (2.3.7). As 
indicated by equation (lw2.5), the coefficient of l/£, in the asymp- 
totic expansion of equation (2.3.8), should be zero. This permits 
eli min ation of one of the- A n in equation (2.3.8). Recall that 
Yi and y^ are boundary edges, that is, points at which the boundary 

conditions change from a specification of chr/dx to a specification 
of cto/Bx. If a boundary edge does not vary with x, the correspond- 
ing singularity in the denomi n ator and the polynomial in the numerator 
are each of one order lower than indicated in equation (2.3.8). The 
symbol y^" again represents points at which there is a discontinuity 

in the v or w boundary condition. If y”' does not vary with x, 

however, the corresponding term does not appear in equation (2.3.8). 
That this is the case may be verified by replacing the discontinuity 
by a continuous change over a narrow interval and then letting the 
interval go to zero. This process does not alter the form of the 
generating function dv/dx and hence this generating function is 
insensitive to discontinuities in the v or w boundary conditions 
provided these occur along lines of constant y. The discontinuities 
are reintroduced into the flow field by the integration with respect 
to x. 


For the wing plan form shown in figure 7(a), with the general 
wing .boundary condition w = f(y) + g(x), equation (2.3.8) becomes 
(for x > c) 


dV 

5E = 


i 


A p + A X E + + jfM £* 

<\A 2 - y 2 2 (E 2 - y 2 2 ) 3 / 2 


i d g( x ) 

dx 


(2.3.9) 


where r = 0 is assumed (which makes A3 = 0) and the Kutta condition 
is imposed at the trailing edge. Equation (2.3.9) is used in sections 
4 to 6 to obtain the aerodynamics of swept wings. 

For a semi-infinite swept wing (fig. 8) with the general wing 
boundary condition w = f(y) + g(x), the generating function is 
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dV 

cix = 


i 


- yi (5 - y2) 3 / 2 _ 


i d s( x ) 
lx 


(2.3.10) 


assuming the Kutta condition at y-^. Moreover, Aj_ = - ^ — (y^+ 3yg) 

so as to satisfy the condition of no net circulation. There is a dis- 
continuity in v across the x-axis (edge of shed vortex sheet) which 
does not appear in equation (2.3.10). This is due to the fact that the 
boundary conditions on both sides of the x-axis are in terms of v 
(i.e., chr/dx = 0) so that the x-axis corresponds to a y“* type 

point which is independent of x. Equation (2.3.10) iB used in 
section 11 to solve unsteady two-dimensional airfoil problems. 


SECTION 3 - INTEGRAL EXPRESSIONS FOR FLOW FIELD 

As mentioned previously, the generating functions Sv/cbc and 
Sv/9^ are used to solve the direct and inverse problems, respectively, 
of slender wing theory. The final solution of a given problem involves 
the integration of the generating function and the elimination of the 
arbitrary A^ so as to satisfy the boundary conditions. The integral 

expressions for the direct and inverse problem are indicated in sec- 
tions 3.1 and 3.2, respectively. The solutions should be considered 
valid for only the upper half -plane (z > 0) because of the arbitrary 
-procedure used to evaluate the branch pointB (section 2.l). The 
solution for the lower half -plane is found from symmetry considerations. 

The boundary conditions to be satisfied are that there be no lift 
acting across the shed vortex sheet and that the specified w distri- 
bution exist over the wing plan form. The condition of zero loading 
across the shed vortex sheet can usually be imposed without difficulty. 
Howevef, in the case of a direct problem, satisfying the specified w 
boundary condition on the wing requires the solution of an integral 
equation. Tiro forms of the integral equation are indicated in 
section 3.3. 

5.1 Integral expressions (direct problem) . - Replacing x by £ 
as an integration variable yields 






W = 


(3.1.1a) 
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(3.1.1b) 


Note that W and V must be continuous functions of x for all three- 
dimensional flow problems where slender body theory Is applicable. 


The complex velocity V is obtained by differentiating W with 
respect to £. A corresponding complex function U can be defined 
as the derivative of W with respect to x. That is. 


U = 


Sw 


= u + 



(3.1.2) 


so that the real part of U is proportional to the loading in the 
z = 0 plane. Moreover, d^W/dx = is = 80 that 



(3.1.3) 


Equation (3.1.3) can be integrated directly and is used to satisfy the 
condition of zero loading in the shed vortex sheet. 


5.2 Integral expressions (inverse problem) . - The corresponding 
expressions for the inverse problem are 



(3.2.1) 
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5.5 Integral equations (direct problem) . - Assume that the cross- 
flow at some section x = c is known. This is true, for example, for 
the swept wing indicated in figure 7. Then, from equation (3.1.1b) 


w(x,y,0) - w(c,y,0) 



(3.3.1) 


where the notation w(x,y,z) means that w is evaluated at (x,y,z). 
If the value of y is such that (x,y,0) corresponds to a point on 
the wing surface (as in fig. 9(a)), then the left-hand side of equa- 
tion (3.3.l) is known and equation (3.3.1) is an integral equation 
for Sv/d£. 


An alternate form of the integral equation can be found by inte- 
grating both sides of equation (3.3.1) in respect to y. Let T) be 
the integration variable in the y-direction. Integrating both sides 
between the limits y and y 2 (x) yields 



where the order of integration on the right-hand side has been re- 
versed. The area of integration is indicated in figure 9(b). Equa- 
tion (3.3.2) can always be reduced to a function of x plus a func- 
tion of y equal to a function of x plus a function of y. The 
functions of x or of y can then be equated, providing an alter- 
nate form of the integral equation for dv/d£. 


Numerical methods are usually required to solve equation (3.3.1) 
or (3.3.2). For the case of a swept wing, the numerical solution of 
equation (3.3.1) requires an integration by parts in order to reduce 
the order of the singularity of the integrand and this introduces 
derivatives of the unknown functions A n into the integrand. How- 
ever, the numerical solution of equation (3.3.2) does not require 
the integration by parts since the order of the singularity # in the 
integrand is reduced by the tj - integration. Thus the numerical 
solution of equation (3.3.2) is usually more straightforward than 
that of equation (3.3.l). The alternate forms of the integral equa- 
tion, for the swept-wing lift problem, are derived in section 4.2. 
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SECTION 4 - LIFT CEF SWEPT WINGS (DIRECT PROBLEM) 

The solution for a swept wing at angle of attack is derived in 
terms of S, a function of x, which must he evaluated from an inte- 
gral equation. Some limiting solutions of the integral equation are 
obtained. 


4.1 Load distribution . - For the lift case, the w boundary con- 
dition on the wing of figure 7 is w = - aUQ. The solution for x < c 
is well known and equals 


W = - iaU Q 



0 


(4.1.1a) 



(4.1.1b) 



(4.1.1c) 


dL j, 
to = 4rt<3ay 2 



(4.1. Id) 


Equations (4.1.1a) and (4.1.1b) are valid at x = c since W and V 
are continuous functions of x. However, dv/dx and dl/to may be 
discontinuous at x = c. 


For x > c, the generating function is (from eq. (2.3.9)) 


EE = ' iaU °l 




+ a 2 £‘ 


aA 2 - n 2 (E 2 - y 2 2 ) 3 /2 


> 


(4.1.2) 


where the symmetry with respect to y has been imposed (i.e., Ap = 0) . 
The functions Aq and Ag are eliminated by satisfying the w bound- 
ary condition on the wing and the boundary condition that there be no 
loading in the wake. 
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The load distribution in the z = 0 plane is proportional to u 
(eq. (l.l.2))3. The value of u on the upper surface of the right- 
hand wing panel is, from equation (3.1.3), 


u = R.P. 



(Ap + A 2 £ 2 )d£ 

(P - y 2 2 ) 3 / 2 


(4.1.3) 


The integration is conducted along the upper surface of the z = 0 
plane with a suitable indentation at yg. If tj is used as the 

integration variable in the y-direction, equation (4.1.3) can be 
written as 


f 

U = <xUq -j- 


(A o + A 2 il 2 )dTi 


Jy 2 


4 


n 2 - y x 2 (y. 


>2 _ ,,2)3/2 


(4.1.4) 


whereby- indicates the infinite part of the improper integral (appendix 
B). Integration of equation (4.1.4) yields (see appendix C) 


u 

-aU 0 - 


Ap + A 2 yx 2 

y 2 (y 2 2 - yi 2 ) 


L 


F(PSk‘) - 


A o + A 2 y 2 2 


y 2 (y 2 - yi ) 


- ^E(p<, 


k’) + 



(4.1.5) 


^Strictly speaking, equation (1.1.2) is applicable when the bound- 
ary conditions are satisfied on the wing surface. When the boundary 
conditions ere satisfied in the z = 0 plane, as is done herein, the 
appropriate expression for pressure is 

( v 2 + w 2 \ 

p - Po = - Po \ nJ3 0 + awU 0 + § ) 

for a configuration at angle of attack a. At any rate, the loading is 
proportional to u, for a wing of zero thickness, since all the other 
terms are symmetric with respect to z. 
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w 




o 


where F(p*,k‘) and E(p*,k') are incomplete elliptic integrals of the 
first and second kind, respectively, with amplitude p 1 and modulus 
k‘ . These have the values 



(4.1.6) 


To satisfy the condition of zero loading in the wake, equation (4.1.5) 
is set equal to zero for y = yq. This gives a relation "between A q 

and Ag. If S, a function of x, is introduced according to 



A o + k &z 
= y 2 (y 2 2 - yi 2 ) 


then equation (4.1.5) becomes 


u = aUo 


(g) 


B(p*,k«) 


F(p l ,k‘) 


El . _y 
K* Y 2 



(4.1.7) 


where K* and E* are complete elliptic integrals of first and second 
kind with modulus k' . Equation (4.1.7) describes the spanwise vari- 
ation of loading on the wing. The unknown function S appears only 
as a scale factor. The generating function can now be written as 


3V , TT (Wl 
35E = " 672 WJ 


f C 2 - Yi 2 

- Y2 2 


E'\ 1 

K 7 V(? 2 - y 2 2 )(C 2 - yi 2 ) 


(4.1.8) 


From equations (1.3.4) and (1.2.5) and the asymptotic form of equation 
(4.1.8), the lift per unit x is, in terms of S, 


dL 

dx 


4jtoqyg (— ) S (l - |i) 


(4.1.9) 


From equations (4.1.7) and (l.3.1l) the suction force at the leading 
edge is 



= - m 2 qS 2 



(4.1.10) 


The function S remains to be determined. 
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4.2 Integral equations for S . - The function S is determined 
from equations (3.3.1) or (3.3.2). Each of these equations is con- 
sidered separately. 

(a) Equation (3.3.l) : Substitution into equation (3.3.1) yields 



where b is the value of yg at x = c and y > b is assumed. The 

path of integration for equations (4.2.1) is indicated in figure 9(a). 
Transforming from £ to yg as the variable of integration giveB 



(4.2.2) 


From equation (4.2.2) it is found that S -*■ 1 for yg -► b. Thus dv/dx 

and dli/dx are continuous functions of x at x = c. It is rather 
remarkable that dL/dx is continuous at x = c since it can be shown 
that there is a pressure discontinuity at this section. Equation 
(4.2.2) is a Volterra type integral equation for S, which must, in 
general, be solved numerically. To permit numerical solution, the 
finite-part operation must be eliminated. This can be done by inte- 
grating the 3/2-order singularity by parts j the result is 
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nr 


o = 


^ 2^2 


.2 _ 


yi 2 dS 


A^/Cy 2 - yi 2 )(y 2 - y2 2 ) 


Sy 2 dy 2 


Ub 


yi 

y2 ^ 2 ; 

(4.2.3) 


It Is noted that the integrand of equation (4.2.3) contains both the 
unknown function S and its derivative. 


(h) Equation (5.5.2) : Substituting into equation (3.3.2) and 

noting that y > b yield 

- Vy2 2 W - b2 + Vy 2 - " b2 



It is desirable to express the right-hand side of equation (4.2.4) as 
a function of x plus a function of y so that the functions of x 
or y can he equated. Taking the imaginary part of the right-hand 
side and utilizing the finite -part technique result in an area of 
integration as indicated in figure 10(a) . This integration can he 
decomposed into two separate integrations as indicated in figures 
10(h) and 10(c). The right-hand side of equation (4.2.4) then 
becomes 
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The first term in equation (4.2.5) is a function of x only while the 
second term is a function of y only. Equating the functions of y 
gives 



Integrating the inner integral and converting from g to yg as the 
integration variable give 



where 

- i) [K - F(p,kj] + E - BO,k) 

( E I \ 

= ggT + ( 1 - gyl F(p,k) - E(p,k) (by Legendre's relation) 

P = jz/j k = yjjyz 

and K and E are complete elliptic integrals with modulus k. Equa- 
tion (4.2.7) provides an alternate Volterra type integral equation for 
S which does not contain the derivative of S in the integrand. 

Equation (4.2.7) differs from equation (15) in reference 11. The 
discrepancy arises from Robinson's treatment of his equation (12). He 
properly states that a function of x can be added to his equation 
(12 ) , but takes this function to be zero. This function of x actu- 
ally does not equal zero except for x = c. Hence Robinson's integral 
equation is correct only for the limiting case of a wing whose trailing 
edge is only slightly swept. This limitation of Robinson’s work was 
previously pointed out by this author in reference 12. 
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4.5 Limiting solutions of Integral equations. - 


(a) Special solution : It has already been established that S = 1 

at x = c. A special solution of equation (4.2.3) is to assume that 
S = 1 for all x. Then 


II - v ^ 

K’ * dy 2 


(4.3.1) 


Integrating 


*1 = k 
b E* - k 2 K* 


(4.3.2) 


which gives the equation of the trailing edge. This result was obtained 
in reference 6, by other methods, and corresponds to the case where 
there is no shed vortex sheet. A more detailed discussion will be 
given in section 7.2. 


(b) k « 1 ; For k « 1, the solution for S is (eq. (D4b) of 
appendix D) 


S = 1 + 4 





dt + ik 2 + 


y In t 2 


0(i4) + 0 t ^l) + 0 © 


(4.3.3) 


where y = dy-^/dyg is the ratio of the trailing-edge slope to the 

leading-edge slope. For y equal to a constant, the integral in 
equation (4.3.3) becomes the logarithmic integral which is tabulated 
in reference 14. 


(c) y : When the ratio of the trailing-edge to leading-edge 

slope is large, the solution for S is (eq. (D5)) 


S = -- -T. 1 (4.3.4) 

a/i - k 2 

which is valid for all values of k. Equation (4.3.4) is equivalent to 
that obtained by Robinson since hiB integral equation is correct when 

X -* m • 


(d) y = 1, k -»• 1 ; This case corresponds to a two-dimensional 
swept wing. The solution for S is S = 1, 
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4.4 Numerical solutions. - The y-coordinate will be nondimension - 
alized with respect to b by means of the notation Y = y/b , Yg = yg/b, 
and so forth. Equation (4.2.7) becomes 




(4.4.1) 


Divide the integration interval into m parts (fig. ll). Let rjg^ n 
represent a mean value between Yg^ n _ 2 _ and Yg^ n . Let T)-j^ n repre- 
sent a mean value between Yp^ n _i and Y]_^ n . When the mean value 
theorem is used, equation (4.4.l) becomes 


Vy 2 - 1 = Sn*n(*2,n - *2,n-l) + 

n=l 

S n ~ (Vy 2 - Y| jn .! - Vy 2 - Y^ n ) (4.4.2) 

where Sn and 3> n are evaluated at i)g^ n . Thus 

P n = ^ 2 ,n/ Y 

^n = ^l^n/ 1 ! 2 ,n 

and so forth. The values of Sj_, Sg, S 3 , . . . are found by succes- 
sively letting m = 1, 2, 3, . . . in equation (4.4-.2). The explicit 
expression for S m is 


m-l 


8 * 


a/y 2 - 1 - T 2,n-l) + ^ y " ^ (a/y 2 - Y§^ - aJi 2 - Y z , n ) 

n=l 

*n(r - l2,B-l) + 1 ni,m a/y 2 - Yl,n-1 


( 4 . 4 . 3 ) 


An alternate expression for Sj^ can be obtained from equation (4.2.3). 

Equation (4.4.3) was evaluated for the 1=1 case with intervals 
(Yg^ n - Yg^ n _ 3 _) = 0.20. The intermediate points T)^ n and T)g, n were 
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taken to be in the middle of each interval. The results are compared 
with equation (4.3.3) in figure 12. The two agree within 2 percent for 
y 2 /b < 2.5. According to figure 12, S decreases from. 1 at y g/b = 1 

to a minimum value of 0.94 at y 2 /b = 2 and then increases slowly to 
its asymptotic value of 1. 

The aerodynamic forces can be found from a numerical, integration 
with the previously found values of S. Thus, the total lift is given 
by 


2itqab^ 


= 1 + 2 



,b 0 A 


S(1 - E*/K' }YgdIg 


m 


= 1 +X S n (l - = 0 /^)^ l,n - Y l,n-l) ( 4 - 4 - 5 ) 


n=l 


bo/b - 1 

where m = r- when equal intervals are used. Similar ex- 

x 2,n " x 2,n-l 

pressions can be deduced for pitching moment, center of pressure, and 
induced drag. Some of these coefficients are presented in figure 13 
for the case y = 1 and dy 2 /dx = constant. 


SECTION 5 - ROLLING SWEPT WING (DIRECT PROBLEM) 

The solution for a rolling swept wing is presented using the pro- 
cedure of section 4. 

5.1 Load distribution and rolling moment . - The boundary condition 
on the wing is w = - (o^y where o>x is the angular velocity. The 
solution for x < c is 


W = - 



V 


% / 2£ 2 - VZ 


i vr- 


2 '4^7 



(5.1.1a) 


(5.1.1b) 
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dV . “x 

S = “ 1 T 


3 r dy 2 

. 72 

)1 

_(t 2 - y 2 2 ) 

372 


dx 


= itq 


“x sf^) 

Uo y 2 3 Ur-J 


(5.1.1c) 
(5.1. Id) 


where cLMx/dx is the ro d lin g moment per unit x. Again., equations 
(5.1.1a) and (5.1.1b) are valid at x = c. 


Prom equation (2.3.9), the generating function for x > c is 

A l£ 


SZ = “ 1 T 


^ 2 - Yj 2 (i 2 - y 2 2 ) 3 / 2 

For a point on the right-hand wing panel 


(5.1.2) 



(5.1.3) 


Define R a function of x, such that 

y2 

Then equation (5.1.3) becomes 


(») 


R = 


y 2 2 - yi 2 


u = T 72 


teW 


y 2 - yi 2 


fy2^ - y 2 

The generating function, in terms of R, is 

i 


(5.1.4) 


He = - 1 ir y 2^2 2 - yi 2 ) R 


a 2 - y 2 2 ) 3 / 2 J 


fe) 


(5.1.5) 
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The rolling moment is (from eq. (1.3.8)). 

ar = ”9 -gj- (y 2 2 - E 

5.2 integral equation for R . - Substituting into equation (3.3.2) 
and equating functions of y yield (for y > b) 





n dii 

Ajn 2 - yi 2 (n 2 - y 2 2 ) 3 / 2 


( 5 . 2 . 1 ) 


Integrating the inner integral in respect to tj, and transforming from 
£ to yg as the integration variable in the outer integral, gives 


io 

p 


& 


y 


'Vjy 2 - b2 


R 




y 2 - yi 2 
y 2 - y2 2 




(5.2.2) 


Fran. equation (5.2.2) it is seen that R = 1 at y 2 = b so that 
dMx/dx and Sv/Sx are continuous at x = c. 

5.5 Limiting solutions of integral equation . - 

(a) k « 1 : For k « 1, the solution for R is (eq. (D7b) of 

appendix D) ~ 

R = 1 + | k 2 + £o(k 4 ) + 0(k 3 /r)] (5.3.1) 

(b) x -» m : In this case, from equation (D8), 

R = — ===== (5.3.2) 

/ \jl - k 2 


valid for all k. 

(c) T = 1, k-> 1 : The solution is R = 2. 

5.4 Numerical solution . - When equation (5.2.2) is nondimension- 
alized with respect to b, it becomes 
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(5.4.1) 


Dividing the integration interval into m parts and using the mean 
value theorem as in section 4.4 give 



m 


- 1 = 






Em = 


- jr B^ 2 - n f jn (a/i 2 - 4,n-l - V? - l|,n) 
n=l 


•^Y 2 - '\|y 2 ~- Y 2 


2,m-l 


(5.4.2) 


Values of R were computed from equation (5.4.2) for the y = 1 case, 
with intervals (Yg^ n - Yg n _q) = 0.20. The results are compared with 

equation (5.3.1) in figure 14. The two agree within 2 percent for 
y 2 /b < 3.0. From figure 14 it is seen that, for y = 1, R increases 

monotonically from a value of 1 at yg/b = 1 to its asymptotic value 
of 2 at yg/b -*■ «*. 


The total rolling moment is 


„ ^ b3 

^ " 1<q "v; x 


b 3 

“’ToT 



H(1 - Yi 2 /Y 2 2 )Y 2 3 “2 


n=l ' J 


(5.4.3) 


The rolling moment for the case y = 1 and dyg/dx = constant 
computed and the results are in figure 15. 


o 

oa 


H 

to 


was 
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SECTION 6 - PITCHING SWEPT WING (DIRECT PROBLEM) 


The solution for a swept wing, pitching about the y-axis, is pre- 
sented following the procedure of section 4. 

6.1 Load distribution . - The boundary condition on the wing is 
w = - ODyX. The solution for x < c is 


V = - 


ioy 


S7 , 

EE = ' ^ 


dL 

dx 


( VF : 

y 2 2 - 



U 2 - 

y2 2 

(&)■ 

■ y2 2 

(£ 2 - y 2 2 ) 3 / 2 

oDy d(y 2 2 x) 


+ ^ 


- 1 


* U 0 dx 

where equations (6.1.1a) and (6.1.1b) are valid at x = c. 

The generating function, for x > c, is (frcm eq. (2.3.9)) 

A 0 + a 2^ 2 + 


(6.1.1a) 

(6.1.1b) 

(6.1.1c) 
(6.1. Id) 


= - iXD y 


- 1 


.V ? 2 - yi 2 (Z 2 - y 2 2 ) 3 / 2 J 


( 6 . 1 . 2 ) 


From equation (3.1.3), the expression for u on the right-hand wing 
panel is 

„ m 1 "*° 1 + yi2yzZ Wn- irO A ° + + yg4(l * k ' 2) pfn- vO Ap + A 2 y g 2 + y 2 * y /y2 - y^ ~] 

L - y i 2 ^ ’ y2(ys 2 - yi 2 ) ’ y2(ys 2 - yi 2 ) y 2 Vy2 2 - y 2 J 


( 6 . 1 . 5 ) 


To satisfy the condition of zero loading in the wake. 


A Q + A 2Nl 2 + yi 2 N2 2 A 0 + A 2y2 2 + N2 4 ( 1 + k ' 2 ) 
^2(y2 2 - yi 2 ) y2(y2 2 - yi 2 ) 


Define a new variable Q such that 
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2y 2 Q 


Aq + A 2 y 2 2 + y^C 1 + k * 2 ) 

y 2 (y 2 2 - yi 2 ) 


Then 


Aq = 2y 2 2 Q(y 2 2 E‘/K‘ - yi 2 ) + y^ 2 
A 2 = 2y 2 2 [Q(l - E-/K«) - l] 


and equation (6.1.3) becomes 


u = 2y 2 o> y Q E(p*,k') - ~ F(P',k') + ^ 


y^ - y^ 


/y g - yi g 

fy 2 2 - y 2 

(6.1.4) 


The generating function, in terms of Q, is 



and the lift per unit x, from equations (1.2.5), (1.3.4), and (6.1.5) 
is 

H = 2nq ^ [4y 2 2 Q(l - E«/K') - (y 2 2 - yi 2 )] (6.1.6) 


6.2 Integral equation for Q . - Substituting equation (6.1.5) into 
equation (3.3.2), utilizing equation (6.1.1a), and equating functions 
of y yield, for y > b. 
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^Vy 2 - 


r w „ r 

v -« Jr 2 fe) 


2y 2 2 Q 


f 1)2 ~ y i 2 eA 

yi 2 - y 2 2 ~ K 7 + 


(*i 2 - y 2 2 ) - y 2 : 


^/(n 2 - yi 2 )(n 2 - y 2 2 ) 


2 fc 2 - n 2 \ 

\ n 2 - y 2 2 / 


an 


( 6 . 2 . 1 ) 


Integrating the inn er integral and transforming from £ to yg as the 
integration variable of the outer integral give 

I " b2 + wfc 1^2 ! ya 2 " 2 jf TayW {^ + T ^ T^a]^ 2 (6 ’ 2 - 2) 

where $ is defined by equation (4.2.7). To investigate equation 
(6.2.2), it is convenient to consider Q as the sum of two functions 

qW and q( 2) defined according to the relation 




(6.2.2a) 


such that 


4 


- b 2 = 



Q 


(1) 


$ + 


^2 / 

ri 


y 2 - yi 2 
y 2 - y2 2 


|dy 2 


(6.2.3) 



dy2 




y 2 - y x 2 


Vy2 _ y 2 2 



q(2) 1$ + 


y 2 

■5‘Y 


y2 _ y, 2 

V-% 1^2 

y - y2 , 

(6.2.4) 


Comparing equation (6.2.3) with equation (4.2.7) shows that qC 1 ) = S. 
Thus only equation (6.2.4) requires further study. It can be shown 
that 


Q = 


1 c f 

2 b \dx ) 


1 + 


cCdyg/dx) 


r 


at x = c so that dl/dx and Sv/Sx are continuous at that section. 

6.5 Limiting solutions of integral equation . - It will be assumed 
that dyg/dg is constant so that dyg/d£ = b/c. 
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(a) k « 1 : From equation (D9b) of appendix D, 

Q = i + [o(k 2 ) + o(k/r)j 

(b) y -> ** : The solution for Q is (eq. (DIO)) 



(6.3.1) 


(•6.3.2) 


which is valid for al 1 

(c) r = 1, k -» 1: 


k. 

In this case, the solution for 


Q is 


Q = 1. 


6.4 Numerical solution . - In nondimensional coordinates, equation 
(6.2.4) becomes 



dY 2 

(dy 2 /dS) 




Y 2 _ Yl 2 


Ji 


Y 2 - Y 2 2 


dY 2 


(6.4.1) 


In some cases, the integral on the left-hand side of equation (6.4.1) 
can be evaluated analytically but it is convenient, and also consistent, 
to evaluate it numerically. From the mean value theorem 



n=l 


(6.4.2) 

0 

The solution for Qj^ 2 ) is then 


o 

3 


to 




OZZ£ 



Similarly., the net lift Is 


^o/ b 


CDyC c. 

2 «* w b 


= 1 - 


(y 2 2 - Y^Jar + \ 


2 ^ ( X - E^)(4,n - t 6 - 4 - 4 ) 


The pitching moment about the y-axis 1 b, for constant, 

p^oA m ^ j \ 

^ - 1 - |f (^2 2 - *1 V “2 + ! 2 ^ ( x - 


3 ^V 0 .g 

2 ** W * ° 


Numerical values of Q were confuted, for the y = 1 case, with intervals (Y 2jn “ ^2,n-l) 

= 0.20. The resultB are indicated in figure 16. According to these results, Q starts fr cm a 
value of 1 at y g/h = 1, increases very slightly and then decreases to a value of about 0.90 

for 3 < yg/b < 4. The curve appears to have reached its minimum and will then presumably in- 
crease to its asymptotic value of 1. 

Numerical calculations of the net lift, pitching moment, and center of pressure are pre- 
sented in figure 17 for y = 1 and dyg/dx constant. 


oj 
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SECTION 7 - LIFT OF SWEPT WINGS (INVERSE PROBLEM) 

In the inverse problem, the shed vortex sheet is specified and the 
corresponding trailing edge is determ i ned. Since the shed vorticity is 
antisymmetric about the y = 0 plane (for a lifting wing) the general 
form for v on the top surface of the shed vortex sheet is 



n=l 


where vp, V2, • • • are constants. The solution corresponding to the 
first term of the exp ans ion is presented in the following sections. 


7.1 Determination of crossflow . - The boundary conditions on the 
upper surface of the z = 0 plane are: 


For 


y x 2 < y2 < y 2 2 

w = - aU 0 


O 

ii 

0 < y2 < y-^2 

v = v i |yj h 


o 

II 

61& 


The boundary conditions are homogeneous for the generating function 
dv/d£« There is a discontinuity in v at y = 0 and the Kutta 
condition is applied at y = iy-^. The generating function is then, 

from equation (2.3.6), 


bV 

5T = 


^0 + J A p£ 2 _ 

0 5 V? 2 - n 2 (S 2 - 


(7.1.1) 


.-and the complex velocity equals 



(Aq + A 2 £^)<1£ 

£ V ? 2 - yi 2 (£ 2 - y 2 2 )^‘ : 


V = - icdJo 


(7.1.2) 
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where Ao and A 2 are functions of x. Equation (7.1.2) can be inte- 
grated directly to give 


V = 


iaU 0 

~2~ 


2(Ap + A 2 y 2 2 ) / k 2 - yp 2 

y 2 2 (y2 2 - yi 2 ) \ v£ 2 - yz 2 



■kp 

yiy2 3 


In 


- ( y^2 + V(e 2 - yi 2 )(C 2 - y 2 2 ) 


r(y 2 - yi)' 


(7.1.3) 


When the appropriate branches are taken, equation (7.1.3) gives V at 
all points in the flow field. The functions Ao and A 2 are elimi- 
nated so as to satisfy w = - aUo on the wing and v = vi|y|/y in 
the vortex sheet. The result is 


A o - - f 5^ w 


A 2 - - (y 2 2 - si ) 1 - i 


1 vi y 2 + yi 2 ^! yi y 2 


dUr 


In 


J 0 y2 - yi * aU 0 y 2 2 - yr 
Substituting into equation (7.1.3) gives 

V 1 


V = - iaU 0 


H 


aU, 


In 


y2 + yi \ £ 2 - yi 2 


y 2 - yi J\\Jz 2 - y 2 2 

1 V 1 , - (yi72 + A Jit 2 - yi 2 )U 2 - y 2 2 )) 2 

e 2 (y 2 - yi) 2 


II - 


jt aU, 


For a point on the wing panel: 

V 1 


V = a,U 0 


m L 2 ~ J^ - 1) 

V* aU 0 y 2 " yi ) 




y 2 - yp 2 


y 2 2 - y 2 


1 y 2 (y 2 2 + yi 2 ) - 2 yi 2 y 2 2 1 

* ^0 COS " y2(y 2 2 _ y^) + *] 


(7.1.4) 


(7.1.5) 
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For a point on the vortex sheet: 

V = v-j_ + IoUq ^ 

y 2 + yi 


i _ji 

re aUn 


yi 2 - y 2 
[ y 2 2 - y 2 


In 


y 2 - yi 


In 


(yiy 2 - Vyi 2 - y 2 ^2 2 ~ y^-y 4 


yi 2 - y 2 
y 2 2 - y 2 


y^(y2 - yi ) 2 


(7.1.6) 


The equation of the trailing edge must now he determined. 

7.2 Equation of trailing edge . - The equation of the trailing edge 
can he found hy calculating the potential at the trailing edge hy two 
different methods and then matching the solutions. The upper surface 
of the z = 0 plane is- considered. 

The potential along the x-axis, for x 5- c, is constant and equals 
oUq! . Integrating in the positive y-direction gives the following 
formula for the potential at the right-hand trailing edge 


<Px = aU o b + Tiyi 


(7.2.1) 


The expression for <p^, obtained hy integrating v from the leading to 
the trailing edge, is 



Equating equations (7.2.l) and (7.2.2) gives 



k 


i ja 
n oXSq 


In 


1 + k 
1 - k 



I 



— — ~ kK'^ 

jt aUo / 


yi 

h 


(7.2.3) 
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Ol 

H 

CO 

O 


where k = y^/yg as before. The right-hand side of equation (7.2.3) is 
a function of k with vq/aUQ appearing as a parameter. If vj/clUq 


is specif ied, equation (7.2.3) defines the trailing edge, 
equation (7,2.3) becomes 


For 


vi 

aU 0 


- 0, 


n = k 

b E* - k^K' 


(7.2.4) 


which was the special result indicated by equation (4.3.2). For this 
case the leading and trailing edges tend to become parallel as x -*■ 00 . 
The spacing is given by lim (y 2 - yq)/b = 2/it. 

x ->.co 

The wing plan forms corresponding to several values of vi/aJJo are 
indicated in figure 18. For all nonzero values of vq/dUo the trailing 
edge has a cusp at x = c (too anal 1 to appear on the figure) . This 
result indicates that when the trailing edge is specified so as not to 
have a cusp (direct problem) , the sidewash v will be zero along the 
x-axis (for all x) but, for x > c, will probably increase quite rap- 
idly with y. Note that for vq / 0, the apex of the trailing edge is 
a triple point for v. As x -*■ “ (or yj)-*-* 0 ), equation (7.2.3) becomes 


Y 1 rt 

aU 0 ~ , 1 + k 2kK 1 

1 - k E* - k 2 K* 


(7.2.5) 


which is an equation for the limiting value of k in terms of vq/aU 0 - 

Thus, as x -*■ 09 , the trailing edge becomes straight (assuming the lead- 
ing edge is straight) and the ratio of the trailing-edge slope to 
leading-edge slope equals the limiting value of k obtained from equa- 
tion (7.2.5). In effect, the flow field tends to become conical as x 
increases. From figure 18 it appears that taking the shed vortex sheet 
equal to Vq |y| /y is adequate to simulate a swept wing having straight 

leading and trailing edges such that y < 1. 

7.3 Lift . - The asymptotic form of equation (7.1.4) is 
aU 0 

t 2 


(l - 1 

^ « aU 0 


In 


y2 + yiV y 2 2 ~ y i 2 
y 2 - yqA 2 


^ 1 V 1 

I - x y l y 2 + • • • 


V = - i 
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The net lift acting on the airfoil is then 
" (y 2 2 ' yi2) 


where y^_ and y 2 are evaluated at x = Cq. 

SECTION 8 - ROLLING SWEPT WING (INVERSE PROBLEM) 

The shed vortex sheet, for a rolling wing, is symmetric about the 
y = 0 plane and v may be expressed as 


i v i ( 2^2 y2 + yi ] 

* aU o \y 2 2 - yi 2 + 111 y 2 - yi / 


ao 



n=l 


where v-^, v 2 , • - . are constants. In the following it will be assumed 
that v = v^ and the corresponding solution will be obtained. 


8.1 Determination of crossflow . - The boundary conditions are that 
w = - (Djj-y on the wing panels and v = yp in the shed vortex sheet. 

The generating function <dV/B£ has the nonhamogeneous boundary condi- 
tion dw/dy = - cDx on the wing panels and therefore must be expressed 
as (eq. (2.3.7)) 


SV _ ^ j Aojb A2^ 2 + 

^ ^ l V ? 2 - yi 2 - y 2 2 ) 3 / 2 

Since there is no net lift, dv/d£ behaves like l/^ 
that 


1 


( 8 . 1 . 1 ) 


t,-* m BO 


A 2 = - \ (yi 2 + 3y2 2 ) 


The function Aq is eliminated by integrating equation (8.1.1) and 
satisfying the w boundary condition on the wing. The result is 


ao - *4 (3n 2 - y 2 2 ) ♦ & 



2 


2 


3120 



3120 


MCA TN 3105 


45 


The expression for v on the wing can then be shown to equal 


- . ,2 
v = -p— y 2 k'^ 


- [p(p*,k*) - E(p',k*)] | - P(p*,k«) - 


fK 2” 

Jy2 - yi 2 

L E ^!] 

Vy 2 2 - y 2 j 


( 8 . 1 . 2 ) 


The value of v at the trailing edge is 


* ^2*“ 
V 1 4 E 


(8.1.3) 


8.2 Determination of trailing edge . - If Vp is specified, equa- 
tion (8.1.3) provides an expression for the trailing edge. Considering 
the limiting case of x = c, however, shows that Vp must equal ax^o/Z, 

which is the value of v for the basic wing at x = c and y = 0. 

Thus vp is no longer a free parameter. The reason for this is as 

follows. The lift is zero along the x-axis of the basic wing. Hence, • 
no modification of the trailing edge, at x = c, can change the span- 
wise lift distribution on the wing sufficiently to change the value of 
the shed vorticity at x = c. Thus v is continuous along the x-axis, 
at x = c, and must equal o> x b/2. The equation for the trailing edge 
is then 


yi _ 2E k 

b it (k«)2 


( 8 . 2 . 1 ) 


which is plotted in figure 19. As x -*• <*, equation (8.2.1) becomes 


lim 

x-*» 


y 2 - yi 

b 


1 

ir 


Equations (8.2.l) and (8.2.2) were previously obtained in reference 7 
by other methods. 

8.3 Rolling moment . - The asymptotic form of is 


dV _ . ^x(y2 2 - yi 2 ) 
8£ 4 


[jy2 2 + yi 2 - 4 (y2 2 - yi 2 ) §] 
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The net rolling moment is then, utilizing equation (1.3.7), 

— ^ = (y2 i 2 y ~ " ^ [ 7y 2 2 + yi 2 - 4 (y2 2 - yi 2 ) |] C 8 - 3 - 1 ) 

rtq uo 

where y^ and yg are evaluated at x = Cq. 

SECTION 9 - PITCHING SWEPT WING (INVERSE PROBLEM) 

The shed vortex sheet is assumed to be v = V]_|y|/y, as in the lift 
problem, and the resulting solution is obtained. 

9.1 Determination of crossflow . - The crossflow is obtained by 
replacing oUq by aye in equations (7.1.1) to (7.1.6). 

9.2 Equation of train rig edge . - The potential along the x-axis, 
for x > c, is cp = ayeb. The potential at the right-hand trailing 

edge can then be expressed as 


<Pl = (DyCb + v-jy-L 


(9.2.1) 


Integrating v from the leading edge gives 




1 V 1 
It OOyX 


In 


1 + k 
1 - k 



k 2 K') + — — (2K* 

1 Jt OiyX V 


'1 

(9.2.2) 


Equating equations (9.2.1) and (9.2.2) and assuming that the leading 
edge is a straight line yg = bx/c give 


y l , 1 7 1 / 1 + k 2kK« \ 

b 2n ca)y y 1 - k jw _ ]j2g-i J 



which is an expression for the trailing edge in terms of k and the 
parameter vj/oay The plan forms corresponding to different values 

of vp/coiy are plotted in figure 20. As in the lift case, the 
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trailing edge las a slight cusp at x = c (except for = 0) . How- 
ever, for x-*«, the trailing edge approaches the leading edge. When 
Vj_ = 0, equation (9.2.3) becomes 

yi = k 

b -s/e* - k 2 K' 


which is one of the curves included in figure 20. 


9.3 Pitching moment. - The net lift acting on the wing upstream of 
a given section x is obtained by replacing oUq in equation (7.3.1) 
by oyx and equals 


L 

2nq U0 


x(y2 2 



I Vj / 27172 Y2 + yj ] 

II \y 2 2 - yi 2 + y 2 - y i)_ 

(9.3.1) 


The pitching moment is obtained by integration. Thus, the pitching 
moment about the leading edge is 



which can be evaluated by means of equation (9.3.1). 


(9.3.2) 


SECTION 10 - WING-BODY COMBINATIONS 

To solve nonplanar problems, such as the flow about wing-body 
combinations, it is necessary to transform the given problem into 
one with planar boundary conditions in order for the previously 
developed generating functions to be applicable. In the following 
section the Joukowski transformation is discussed. The solution 
for the flow about a highly swept wing mounted on a circular cylin- 
der is then indicated. 

The x,y,z-coordinate system is considered, herein, as a body axis 
system rather than a wind axis system in order that the developments 
parallel those of the isolated swept wing problem. The perturbation 
velocities u, v, and w are considered as parallel to the x, y, 
and z coordinates, respectively. This creates no essential change 
except for the fact that the pressure formula (eq. (1.1.2)) becomes 
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P - Pq PQ 



4- ccwUq + 



when the hody axis is at angle of attack a in respect to the free 
stream. Also, the lift and moment formulas (eqs. (1.3.3) to (1.3.5) 
and (1.3.7)) must he evaluated after transforming from the body axis 
system to a wind axis system. 

10.1 Joukowski transformation . - Consider the problem of a highly 
swept wing on a circular cylinder of radius a (figs. 2l(a) and 21(b)). 
The Joukowski transformation 

T = 5 - a 2 /5 (10.1.1) 


transforms the configuration such that the body becomes a vertical cut 
(of width 4a) on the z-axis of_the ^-plane (fig. 21(c)). The wing 
panels remain as cuts in the z = 0 plane. The velocities at cor- 
responding points in the two planes are related by 


V 



( 10 . 1 . 2 ) 


where V = v 


iw is the complex velocity in the ^-plane . But 


where 


dT 

dT 

= 1 + *! 
5 2 

- 0 *( 

§) 2 c ° e ze ] - 1 0 

sin 20 

and 0 

are defined by the relation 5 = 

re 1 ®. Thus 

V = V 

t * (# 

> - 

cos 26 

- w sin 20 


w = w 

Hs) : 

cos 20 

| + v sin 20 


!t> 

1 l 

11 

£ 

cos 9 + 

w sin 0^ 


-) 2 sin 0 


(10.1.3) 


where v r is the radial velocity in the 5 -plane. Equations (10.1.3) 

are particularly useful for relating the velocity boundary conditions 
in the and ^-planes. 


Since body axes sire being used, the boundELry condition on the body 
surface, in the £ -plane, is that v r = 0. The corresponding boundary 

condition_in the ^"-plane is that v = 0 along the vertical, cut 
(- i2a < £ = iz < i2a). This boundary condition is automati cally 
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M 

O 
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I 


8 


satisfied, by symmetry, for all problems when tbe _w boundary condi- 
tions on the wing panels are symmetric about the y = 0 plane (i.e., 
lift anrl pitch problems) . For these problems, only the boundary con- 
ditions in the z = 0 plane and at infinity require attention. Thus, 
the wing-body problem is transformed, by equation (lO.l.l) into an 
equivalent isolated swept-wing problem (for the symmetric case)^ . 

The asymptotic form of the complex velocity in the ^-plane is 


V = - 


ay 


+ 0 


1 

_«) 3 _ 


The lift in the T^-plane is then (from eqs. (1.3.4) and (1.3.3)) 


dL 

dx 


= i.P.tertpoUQdfy/dx) 


(10.1.4a) 

(10.1.4b) 


Equation (10.1.4b) corresponds to a configuration which is pointed at 
its upstream end. The complex velocity in the physical plane is 



Note that the coefficient of the leading term in the asymptotic expan- 
sion for V is unaffected by a transformation from body axes to wind 
axes. The lift per unit x in the physical plane is then (eq. (l.3.4)) 


dL _ __ 

— = I.P. 

dx 


dfi dM CB 

2 * p 0 U 0 + p oV — ^2 


(10.1.5) 


But £g = - lax + constant and A cs is constant since a cylindrical 
body is assumed. Equation (10.1.5) then becomes 


■^When w is not symmetric about the y = 0 plane (as in a roll 
problem) the boundary conditions along the vertical cut are not auto- 
matically satisfied. In these cases it may be advisable to use the 
transformation 1 = ^ + a2/£ which transforms the circle into a cut 
(- 2a < X = y < 2a) along the y-axis. Thus the problem is transformed 
into one wherein the boundary conditions are specified along the y-axis 
and at infinity. The problem is now completely planar but involves 
three "panelB." In some cases it can be handled by the methods of 
section 2. 
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= I • P • ( 2itp 0 U 0 df 2 /dx ) 


(10.1.6a) 


L = I.P.(2rtp 0 U 0 f 1 ) (10.1.6b) 

Equation (l0.1.6b) represents the net lift of the configuration in the 
physical plane provided the body remains cylindrical upstream of the 
wing. For pointed forebodies, the additional term rt Po^0^ a ^' a '> repre- 
senting forebody lift, must be added to the right-hand side of equa- 
tion (lO.lJSb). Compar ing equations (10.1.4b) and (l0.1.6b) shows 
that L = L. That is, the lifts of the configurations in the physical 
and transformed planes are equal (providing the body remains cylindri- 
cal upstream of the wing) . This result was derived in reference 15 
by another method. 

10.2 Lift of swept wing on cylindrical body . - A swept wing is 
mounted on a circular cylinder and the configuration is at angle of 
attack a. The x,y, z-coordinate system is based on the body axis, 
and the following transformations are made: 

(1) The flow is transformed from the £ -plane to the £ -plane by 
equation (lO.l.l). 

(2) A uniform flow w =±= - aU 0 is added so that the configuration 

is translating, in the £ -plane , with velocity w = - aUQ in a fluid 
otherwise at rest. 


The solution to flow ( 2 ) can be found by the methods of sections 4 and 
7. The solution to the original problem is then found by reversing 
this procedure. 


(a) Direct problem : The generating function for flow ( 2 ) is, frcm 

equation (4.1.8), 


3V , TT - d ^2 Q 
31 = " iaU 0 y 2 dx~ S 


"(I ) 2 

_(D 2 


(D 2 - C^) 2 f. 


- fra) 2 K 'J V[CD 2 - - (Tl) 2 ] 


( 10 . 2 . 1 ) 


where K* nnfl E* are complete elliptic integrals of first and second 

kinds with modulus k* = 1 - (y-^/y g ) 2 . The lift of the configuration 

is (eq. (4.1.9)) 


dL . 

EE = 4mqy ' 


- ( E‘\ 

^2 Vdx / S V 1 " kJ 


( 10 . 2 . 2 ) 
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The unknown function S is found from the integral equation 
(eq. (4.2.7)) 


CM 

H 

CO 

o 


where 


a/j) 2 - (b) 2 = 


p 

_ y 2 /(y ) 2 - (?i ) 2 

S dy 2 
b 

y V(y ) 2 - (y2) 2 _ 


(10.2.3) 


5 =(|r - l)[K - F(p,k)] + E - E(p,k) 

p = Y2/y 
S = Yi/y2 


and K and E are complete elliptic integrals with modulus k. Lim- 
iting solutions for S are given in section 4.3. The solution of the 
original problem is found by reversing steps (2) and (l). Note that 

T = 5 - a 2/5 

b = b - a 2 /b 

7i = yi - a 2 /yi 
y 2 = y 2 - a 2 /y 2 
u = u 

SV SV If df \ 

~Sx = SE/^dt J 

L = L 


The details need not be given. 

(b) Inverse problem : From equation (10.1.2), the relation between 

the shed vortex sheets in the physical and transformed planes is 


v 


v 
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Therefore the sidewash on the top surface of the shed vortex sheet in 
the physical plane will be assumed to be 


V = Vj_ 



!XL 

y 


where vp is a constant. The shed vortex sheet in the transformed 
plane is then 


V = V]_ 


in 

y 


and the solution follows directly frcm the results of section 7. 


SECTION 11 - UNSTEADY TWO-DIMENSIONAL INCOMPRESSIBLE FLOWS 

The flow field due to the motion of a two-dimensional body in an 
incompressible fluid, otherwise at rest, is discussed. The application 
of generating functions for the solution of unsteady airfoil problems 
is indicated. 

11.1 General considerations . - The equations which arise in studies 
of unsteady two-dimensional incompressible flows are closely analogous 
to those employed in slender body theory. The analogy between these two 
classes of flow will be established by comparing the classical equations 
associated with two-dimensional unsteady incompressible flows to the 
equations, derived by Ward, for the flow around a slender body. 

Assume a two-dimensional coordinate system, fixed in space, such 
that the fluid, far from the body, is at rest. The velocity potential 
satisfies 


5 2 <p d 2 (p 

cfy -2 £ z 2 


( 11 . 1 . 1 ) 


where q> contains the time t as a parameter. The velocity potential, 
in this case, is not a perturbation potential and the velocities v 
and w represent the net velocities in the flow field. The boundary 
condition on the body is that 


v n = 


dv 

dt 


( 11 . 1 . 2 ) 


where v n is the velocity normal to the body surface and v is the 
normal coordinate of an orthogonal coordinate system chosen to be 
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normal and tangential to the body surface at a given instant. The pres- 
sure at any point, from the Bernoulli equation for unsteady flow, is 


P - p(t) = - P 0 < 


Sep 

St + “ 


1 r ?^ 2 fa) 2 

2 I \.Sy/ 


(11.1.3) 


where p(t) is an arbitrary function of time. The Laurent expansions 
for W and V are 


W = f 






m 


(11.1.4) 


m=l 


V - S (tT - ir ) r - £ (U.1.5) 

m=l 


where f, f-j_, f 2 , . . . are functions of t. The impulse required to 

generate the motion, at any instant, is found from a contour integral 
about the body and equals (assuming T = 0) 



( 11 . 1 . 6 ) 


Equation (11.1.6) is a classical result derived by Kelvin. Neumark 
(ref. 16) appears to have been the first to evaluate equation (11.1.6) 
by replacing cp by W and using Cauchy's theorem. However, Neumark 
solved the special case of a nondeforming body. For the general case, 
following the procedure used to obtain equation (1.3.3) from equation 
(l.3.l), equation (U.1.6) becomes 


!y + il z = 2np 0 



1 ^(•^•cs^g) 
+ 2n dt 


(11.1.7) 


The force per unit span (ly + i Z z ), at 


any instant. 


is then 


l 


y 




( 11 . 1 . 8 ) 
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Equations (ll.l.l) to (11.1.8) are similar to equations (1.1.4), 
(1.1.5), (1.1.2), (1.2.2), (1.2.3), (1.3.1), (1.3.3), and (1.3.4). 
The two sets of equations become identical if the transformation 
x = Uot + constant is introduced provided the arbitrary functions 
p(t) and f in equations (11.1.3) and (11.1.4) are taken equal to 
p 0 and f in equations (1.1.2) and (1.2.2), respectively. The 

equivalence of equations (11.1.6) and (l.3.l) is established by 


Iy + U z 



[ly + iiz)dt - — 


0 



. . Fy + iFz 

(2y + iZ z )dx = y 


Thus, the steady three-dimensional flow about a slender body can be 
transformed into an equivalent two-dimensional unsteady flow about a 
cylinder whose cross section varies with Uot in the same way that 

the cross section of the original slender body varies with x, and 

vice versa, provided p(t) and f (eqs. (11.1.3) and (11.1.4)), 
which are essentially boundary conditions, are taken equal to Pq 

and f (eqs. (1.1.2) and (1.2.2)) 5 . Note that p(t), p 0 , and f 

do not contribute to the lift. Muhk used this equivalence in his 

studies of airships and many researchers have since referred to it. 

The generating function approach can thus be used for solving 
unsteady two-dimensional airfoil problems. The generating function 
is now Sv/St. When v is used as the integration variable for t, 
equations (3.1.1) became 


^Another approach to this analog is to consider the differential 
equations and boundary conditions for the problem of two-dimensional 
■unsteady motion in a slightly compressible fluid. In this case p(t) 
and f are no longer arbitrary boundary conditions and the analog 
with the three-d im ensional steady flow past a slender body is even 
more striking. For example, the wave drag of the slender body cor- 
responds to the energy radiated as sound in the equivalent unsteady 
two-dimensional problem. 
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Moreover, U is now defined as U = Sw/Bt so that equation (3.1.3) 
becomes 


U 




( 11 . 1 . 10 ) 


One complication arises; namely, a finite impulse can be imparted to the 
system so that ¥ and V may be discontinuous functions of t. Such a 
discontinuity occurs at t = 0 for the problem discussed in section 
11.3. The integrals in equations (11.1.9) must then be considered as 
Stieltges integrals. 


Assuming that the flow at time t = t a is known, the integral 
equations, corresponding to equations (3.3.1) and (3.3.2), are 


dV 


w(t,y,0) - v(t a ,y,0) = I.P. - I ^ 


dr 


z=0 


■J2< x ) 


w(t,Ti,0)dn + [jf(ta,y (t),0) - T|r(t a ,y,0)j 


( 11 . 1 . 11 ) 


= I.P. 



* y 2 (t) 


L. 


a u 


(£) 


dtl 

— >z=0 


( 11 . 1 . 12 ) 


11.2 Two-Dimensional airfoils . - The flow about a zero thickness 
airfoil, moving with velocity Vq along the y-axis in the positive 

y-direction, will be discussed. For convenience, let the chord of 
the airfoil be b and let its trailing edge be the origin of the 
coordinate system at t = 0. If, at t = 0, the airfoil starts to 
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move with velocity Vo, it occupies a strip of width b in the ty-plane 
(fig. 22). The w boundary condition on the strip is defined by the 
prescribed motion of the airfoil. There is a vortex sheet of unknown 
strength behind the trailing edge. Therefore, the generating function 
cW/Bt is applicable. For many problems the w • boundary condition on 
the wing is of the form w = f(y) + g(t). The generating function for 
these cases can be written, from equation (2.3.10), 


dV 

3t = 


i 


An - 


flg(t) 

dt 


(yq + 3y 2 ) + 

VS - yi (S - V2) 5 / 2 


ag(t) 

dt 


dg(t) 

dt 


( 11 . 2 . 1 ) 


where the Kutta condition has been applied at y^ and T = 0 is 
assumed. If w = f (y), the generating function is 


dV ^0 

35 vr^i u - 


( 11 . 2 . 2 ) 


The problem of an airfoil starting impulsively from rest, and maintain- 
ing a constant angle of attack, can be solved by equations (11.2.1) and 
(U.2.2). (if the airfoil is accelerating, dg(t)/dt = - adV 0 /dt in 
eq. (H.2.1), and if the airfoil moves with constant velocity, eq. 
(U.2.2) is used. ) 


11.5 Wagner problem . - The problem of an airfoil, starting impul- 
sively from rest, and moving with constant velocity and angle of attack, 
was first solved by Wagner (ref. 17). Recently, it was discussed, from 
the point of view of slender wing theory, in reference 18. The present 
approach differs from reference 18 in that the generating function is 
used to formulate the problem. 

If the motion is considered to start impulsively, at t = 0, as 
indicated in figure 23, the crossflow is discontinuous across the 
y-axis. Approaching the y-axis frcm t < 0 gives a zero velocity 
field. Approaching the y-axis from t > 0 gives a velocity field 
equivalent to that about a flat plate translating with velocity 
w = - aV 0 . The solution for this flow field is 

W = - iaV 0 O - b) - (£ - b/2)] 

V = - iaVA I £ ~ b / 2 - 1 

0 |_V£(£ - b) 


(11. 3. la) 
(11. 3. lb) 
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i 

B 


For t > 0, tie generating function is 


dV 

= 


- idVp 


^0 

V£ - Yi .(£ - J2) 5 ^ 2 _ 


(11.3.2) 


The loading on the wing is proportional to bty/bt. For a point on the 
upper wing surface, frcm equation (ll.l.lO), 


-aToAo 


■r 






Vn - yi Cy2 - n) 3 / 2 


^cAo p - 

* vy2 - y 


(11.3. 3) 


Substituting equation (11.3.2) into equation (10.1.12), with t a = 0, 
utilizing equation (ll.3.1a), and equating functions of y yield 


Jj2(y) 


Vy(y - b ) = - 



ApdT] 

Vn - yp fa - Y2) 3 / 2 


Integrating the inner integral and transforming from r to y 2 as the 
integration variable yield. 


- b) m f (^) ^2 (11.3.4) 

Equation (11.3.4) is an integral equation for Ap. The lift on the air- 
foil, from equations (U.1.5) and (11.1.8) and the asymptotic form of 
equation (11.3.2), is 

( 2Ap\ 

— j (11.3.5) 


It can be shown that 


2A 0 

bVp 


for 

2* 


t = 0, and increases monotonically 


to the value 1 for t Equation (11.3. 4) has been inverted,, by the 

Laplace transform, and the results are identical to those of 
reference 19. 
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It is noted that the problem is formulated, from the beginning, in 
terms of lift, which is the variable of primary interest. In reference 
18, the shed vortex sheet is first determ in ed and an additional inte- 
gration is required to obtain the lift. 


SECTION 12 - CONCLUDING REMARKS 

The use of the generating functions and ciV/clx has been 

described for a wide variety of applications. 

Problems for which dv/5£ is applicable may possibly be solved, 
with equal facility, by other methods. But, the use of the generating 
function Bv/Sx seems to have several advantages over other possible 
approaches to the direct problem of slender wings having swept trailing 
edges. First, the dependence of the flow on upstream conditions is 
initially removed an d certain general features of the flow field can 
be determined immediately. Thus, for the lifting and rolling wing, 
the pressure distribution at each chordwise station can be expressed, 
to within a scale factor, without consideration of upstream conditions. 
Second, the problem is formulated directly in terms of quantities which 
define the lift and moments. Finally, the method delays until the last 
stages the problem of solving an integral 1 equation. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 14, 1953 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 




functions of x (or t) 

CM 

H 

DO 

A cs 

cross-sectional area 

O 

Ap f 

plan form area 


a 

radius of circular cylinder 


h 

semi span of swept wing at x = c (fig. 7) 


b 0 

maximum semispan of swept wing (fig. 7) 

M 

o 

aJ 

• • • 

constants 

<■£1 

CO 

i 

k. 

C L 

lift coefficient (= L/q Ap f ) 


Cl 

rolling moment coefficient (= Mx/q hoApf) 


Cfa. 

pitching moment (about x = 0) coefficient 

(= V 9 c 0fpf) 


c 

root chord of swept wing (fig. 7) 


c 0 

crver-all length of swept wing (fig. 7) 


D^Di^Dg^ . . • 

constants 


D i 

induced drag 


E 

complete elliptic integral of second kind with 
modulus k 


E* 

complete elliptic integral of second kind with 
modulus k 1 


E(M=) 

incomplete elliptic integral of second kind with 
amplitude p and modulus k 


l*0,k) 

incomplete elliptic integral of first kind with 


amplitude {3 and modulus k 
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I.P. 


K 


K* 


k 

L 

ly } l Z 

Mq 

P 

P0 

Q,Q(l),Q(2) 

q 

R 


R.P. 


S 

t 

U 

U( 

u 


net forces in x,y, z-directions, respectively 

functions of x (or t) 

imaginary part of complex function 

complete elliptic integral of first kind with 
modulus k 

complete elliptic integral of first kind of modulus 
k' 

7±/72 ( k * = V 1 " k2 ) 

net lift force acting on configuration 
sectional force in y, z-directions, respectively 
moment about x,y, and z-axes, respectively 
free stream Mach number 
pressure 

free stream pressure 

functions of x arising in solution for pitching 
swept wings (section 6) 

dynamic pressure (pqUq 2 /2) 

function of x arising in solution for rolling 
swept wing (section 5) 

real part of complex function 

function of x arising in solution for lift of 
swept wing (section 4) 

time 

derivative of W in respect to x (or t) 
free stream velocity 
perturbation velocity in x-direction 
complex velocity (v - iw) 


V 
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v 0 

y 

W 

w 

x 

x c.p. 

y 

z 

a 

0 


0 * 


r 

r 


£ 

x 

? 

Po 

T 

$ 

<P 


flight velocity in y-dlrectlon 
perturbation yelocity in y-direction 
potential function (= <p + ity) 
perturbation yelocity in z-direction 
coordinate axis parallel to free stream 
center of pressure 
coordinate axis in span direction 
coordinate axis 
angle of attack 


y2 /y 



net circulation in yz-plane 

ratio of trailing edge to leading edge slope 
( dyp/dxN 

\ dy 2 /dxj 
y + iz 

centroid of cross-sectional area (= jg + iZg) 

integration variable in y-direction 
In (4/k) 

integration variable in x-direction 
free stream density 
integration variable for t 
see equation (4.2.7) 
perturbation yelocity potential 
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CDXjCDy 

Special notation: 

s 

n 

y = yi(x) 

x = x-L(y) 

y = y2 Cx)l 

= x 2 (y)j 


crossflow stream function 

angular velocities about x~ and y-axes, respectively 

finite part of improper integral (appendix B) 

quantity in plane 
equation of trailing edge 

equation of leading edge 


x 
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APPENDIX B 

FINITE PART OF IMPROPER INTEGRALS 

Many of the integrals occurring in. the body of the report are of 
the form 

(Bl) 

to be evaluated for y < yg. The path indicated in the following fig- 
ure is chosen 


In = R.P. 



f(P _ 
'(5 - y 2 ) 3 /2 




Z* 


(yjQ) « « « 

^2 ^ 


(assuming f(£) analytic along the path), and equation (Bl) becomes 


I]_ = lim 
e-K) 


R.P. 


^ *<» , ,, . r 2 " 

(5 - Y 2) 3 / 2 (t - Jz ) 5 / 2 

vjy 2+e 


d-n 


(B2) 


Letting £ = yg + in the first of the integrals, integrating, and 

taking the real part give 


I 


1 


lim 

e-K) 


2f(y 2 ) 

VF 


2+e 

fM , 
(n - yg ) 3 / 2 



(B3) 


Equations of this type (eq. (Bl)) occur quite frequently and the ex- 
plicit representation of equations (B2) and (B3) becomes tedious. 
Therefore, to economize, the "finite part" concept is introduced. 

The finite part of an improper integral, having a 3/2-order singu- 
larity at a limit of integration, is defined by 
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— f (S) , a.u. 

- 2f(a) 

^a-e 

P f(€) 

(a - 5 )3/2 ^ ^ 

J 

1 (a - S)3/2 


- 2f(a) 

° f(€) 

(5 - a)3/2 45 ' ^ 

VT 

(€ - a) 3 / 2 


(B4) 




(Ja+6 


Equation (Bl) can then be written 

-72 

It = 

Similarly, for y < y 2 . 



(^- C y 2 )3/2 d71 


(B5) 


i.p. ' y — urn , = f _Jhl 




(£ - Y 2) 3 / 2 f ( y 2 - t }) 3 / 2 

^2 


dT] 


(B6) 


Thus, the finite-part technique is essentially mathematical short-, 
hand, in the present report, since it avoids an explicit representation 
of the limiting processes required to obtain the real or imaginary part 
of complex integrals. The finite-part technique can be generalized and 
has many applications in aerodynamics (see, for example, refs. 20 
and 21) . 
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APPENDIX C 


ELLIPTIC INTEGRALS 


(a) Interval ^ y^ " : Let 

0 = ?zh * = y-Jyz 



Then: 
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(b) Interval y^ yg: Let 
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w 

H 

DO 

O 


n Y 2 


n 4(i n 


^2 


•Vn 2 - yi 2 (y 2 2 - i 2 ) 3 / 2 k ’ 2 


k 2 F(p‘,k*) - (1 + k' 2 ) E(0*,k«) - 


/2 _ v-i 2 




*2 V y 2 2 - y 2 


(c) Expansions : For k^ « 1: Let 

X = In (4/k) (k = Yi/y2) 

(p = y 2 /y) 


Then: 


<p = sin"- L 0 


K» = X + X - 1 k 2 + Jj 
4 64 


k 4 + 


E» = 1 + 

El 

K« 


k* + . . . 


K - f H + 2 {4 ) + 9 Or ) 2 + • • ] 

E ' f t - 2 Of) - 3 Of ) 2 + ■ • •] 

(- 1 ) 

( x - 1) k2 + ( x - i) 

= x [ 1 + (i-i + s ) k2 + °( k 4 i] 

K - F(0,k) = K ^1 - <p^ + sin <p cob cp |jE k 2 + 0(k 4 )j 

E - E(p,k) = E ^1 - E <p^ _ ein cp CO s cp |ji k 2 + 0(k 4 )j 
$ = (cos <p)/X + 0(k^ cos cp)' 
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APPENDIX D 


LIMITING SOLUTIONS OF INTEGRAL EQUATIONS 

Same limiting solutions of the integral equations occurring in 
sections 4 to 6 are obtained. The coordinates are nondimensionalized 
with respect to b so that Y = y/b, Y 2 = y 2 /b Y^_ = y q/b. 

(a) Lift case : Equation (4.2.7) may be written 


rv 


-\/y 2 - 1 = 


Ul 


S dYo $ 


• t { 


Y 2 - Yq 2 
Y 2 - Y 2 2 


(Dl) 


For yi /y 2 = k « 1, # = Vy 2 - Y 2 2 /YX + O^k 2 'n/y 2 - Y 2 2 /y), from 
appendix C. Equation (Dl) can then he written as 


0 = 


dYo S 


S -)/Y 2 - Y 2 2 


YX 


SY 2 

— v 


Y c - Y-i 


Y 2 - Y 2 2 ^/ y 2 _ y 2 


Jl 


0 "\/ y2 - y 2 2 ^ 


Integrating the term containing X by parts yields 



72^*2 



XY ^ 2 + 


Vy 2 - y 2 2| 

S r? o f Sk 2 (y 2 - y 2 2 )"\ 

- Vy 2 - Y X 2 - 1 + 0 ^ y ) 


(D2) 


Define r = dyq/dy 2 , where y is the ratio of the trailing-edge slope 
to the leading-edge slope. Then, for k<< 1, 


o 
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Y 2 = 1 + 


. X + 0 (* 

r - k W 


“2 = Pi^ - v C 1 + °( k A-)] 


y - k r 


Assume S = 1 + 0 
(D4) ) . Then 

nY2 


( r ~ ^i k ) + which will he verified later (eq. 

r\k/ 4 


4 dl 2 = ‘ 4 — y~T + 0 ( ~ p ~ k - — ) +0 (-— T~v 

XY 2 T In t y r 2 in k / \T In k 

do 




Sk 2 (Y 2 - Yo 2 ) 




Vi - (yi/y) 2 •= Vi - k 2 + o(k 3 /y) 

Equation (D2) can then he written 


^2^2 


k/4 


- 4 


dt 


VY 2 - Y 2 2 T 11115 


+ s 


VTT7 2 -i + o(^) + o(-^-)l ( D3 ) 


The solution of equation (D3) is obtained by setting the integrand 
equal to zero. The resulting expression for S is 


S = 


•“Vl - k2 


1 + 4 



k/4 


dt 


r la t 


+ 0 


+ 0 ( >e!_) 

\r J \r 2 in k/ 


,k/4 


= 1 + 4 


— ^~r + ^ k 2 + 0(k 4 ) + 0 
r In t 2 v ' 


(D4a) 


©*»(rfr) 

(D4h) 


which is the solution for S valid for small k. 
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For the case where the trailing edge is only slightly swept 
(t -*■ *° ) j equation (D4a) gives 


S = 


"Nil - k 2 

which can he shown to be valid for all k. 

(b) Roll case: Equation (5.2.2) can be written 



1 + 0 


m 


dJ 2 


R = 




(D5) 


(D6) 


A/l - k 2 

= 1 + | k 2 + 0(k 4 ) + 0 
For r-+ m y equation (D7a) becomes 

R = 

'yjl - k 2 

which can be shown to be valid for all k. 

(c) Pitch case : It will be assumed that dy 2 /di is constant. 

Then equation (6.2.2a) becomes 


(D7a) 

(D7b) 

(D8) 


Q = QW + q( 2 ) 


where = S. The integral equation for Q^ 2 ) may be written, from 

equation (6.2.4) and with the expanded form for <i> y 
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the solution for q(2) 

is then q( 2) = i + o(k/r). 
= ■ ^ 1 ^ + 0 k ) j then, for small k. 



But 


= 1 + 0(k 2 ) + 0 ^ 
For x_ -*• » } equation (D9a) yields • 



valid for all values of k. 


(D9a) 

(D9b) 

(DIO) 
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APPENDIX E 


COMPARISON WITH REFERENCE 13 

In reference 13, Mangier finds the lift, roll, and pitch solutions 
for highly swept wings (direct problem) by a method similar to that of 
the present report. The two papers are compared herein. 

From equation (3.1.3) 


3v _ 3U _ 1 3U 
35 i 5z 


(El) 


so that the generating function 3v/3x can also he considered as ~ 

i oz 

1 

Mangier writes, by inspection, the expression for j ^ for the lift. 


roll, and pitch problems. Each is in terms of a single unknown function 
of x which is determined from an integral equation. The integral 
equation is, in effect, the one that is obtained by using 
equation (3.3.1). 


The basic lift solution in both papers is identical. Mangier' s H 
is the same as the S introduced in section 4. The pitch solutions are 
also In agreement. Mangier* s H^ is related to the Q of section 6 by 


Q = 


1 

2 


xHq dy 2 

„ y 2 



(E2) 


His integral equation for the pitch case can be obtained from equation 
(3.3.1) but does not appear in the present paper. For the roll problem. 
Mangier * s generating function is the same as that used herein and his 
Hp is identical to the R of section 5. However, Mangier appears to 

have made an error in determining his final integral equation for the 
roll case. His equation (72) should be (in the present notation) 



3/2 0 
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Equation (E3) differs from Mangier 1 s equation (72) in that the left- han d 
side of the latter is zero. Equation (E3) can he obtained from equation 
(3.3.1) or by the direct integration of Mangier' s equation (7l), using 
R = 1 fo^ -$2 b. (Gibe numberless equation which precedes equation 

(72) in Mangier' s report is valid only for y < b and is thus used 
incorrectly by Mangier to obtain his equation (72).) Equation (E3) Is 
considerably more complicated than equation (5.2.2) in section 5. 

The treatment of the swept wing problem in the present paper dif- 
fers from Mangier* s in several other respects. Eor example, explicit 
expressions for dv/Sx are given herein which permit the solution of 
all problems which can be handled by the generating function approach. 
(Mangier does not discuss how he goes about getting his expression for 
the generating functions. However, he does show that they satisfy the 
boundary conditions and are unique.) Also, the forces and moments axe 
evaluated, herein, by considering the leading terms in the asymptotic 
expansion of the generating function. This is considerably simpler 
than integrating pressures over the wing surface, as is done by Mangier. 
On the other hand, Mangier has presented numerical results for a wider 
range of plan forms. 
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(b) Crossflow. 


Figure 1. - Slender body in free stream U Q . 


3120 



MCA TN 3105 



Figure 2. - Contours for obtaining forces on body 
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(a) Uniform translation of m wing panels without intermediate vortex 
sheets . 



c Hr 

dv 

Sw 


5v 

dw ^ 






0 




v=0 

w^C-^ 

v=Dj 

W = 

C 2 

V=D m-l 

*= C m 

v=0 


'2m 


00 m Wing panels with constant values of v and w specified in 
z=0 plane . 



(c) Discontinuity in w boundary condition on wing panel. 

Figure 6. - Notation and boundary conditions for solution of multiwing 
panel problems . 
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Figure 18. - Trailing-edge shapes (inverse lift 
problem) . 
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Figure 19. - Trailing-edge shape (inverse roll problem). 









